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Abstract 

We describe SPICE: Simulation Package for Including Flavor in Collider Events. SPICE takes as 
input two ingredients: a standard flavor-conserving supersymmetric spectrum and a set of flavor- 
violating slepton mass parameters, both of which are specified at some high "mediation" scale. 
SPICE then combines these two ingredients to form a flavor-violating model, determines the result- 
ing low-energy spectrum and branching ratios, and outputs HERWIG and SUSY LesHouches files, 
which may be used to generate collider events. The flavor-conserving model may be any of the 
standard supersymmetric models, including minimal supergravity, minimal gauge-mediated super- 
symmetry breaking, and anomaly-mediated supersymmetry breaking supplemented by a universal 
scalar mass. The flavor-violating contributions may be specifled in a number of ways, from spec- 
ifying charges of flelds under horizontal symmetries to completely specifying all flavor-violating 
parameters. SPICE is fully documented and publicly available, and is intended to be a user-friendly 
aid in the study of flavor at the Large Hadron Collider and other future colliders. 

PACS numbers: ll.30.Hv, 12.15.Ff, 14.60.Pq, 12.60.Jv, 13.85.-t 
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I. PROGRAM SUMMARY 



Program title: SPICE 

Programming language: C++ 

Computer: Personal computer 

Operating System: Tested on Scientific Linux 4.x 

Keywords: Supersymmetry, flavor physics, lepton flavor violation 

PACS: ll.30.Hv, 12.15.Ff, 14.60.Pq, 12.60.Jv, 13.85.-t 

External routines: SOFTSUSY Q and SUSYHIT Q 

Nature of problem: Simulation programs are required to compare theoretical models in 
particle physics with present and future data at particle colliders. SPICE determines the 
masses and decay branching ratios of supersymmetric particles in theories with lepton 
flavor violation. The inputs are the parameters of any of several standard flavor- conserving 
supersymmetric models, supplemented by flavor-violating parameters determined, for 
example, by horizontal flavor symmetries. The output are flies that may be used for 
detailed simulation of supersymmetric events at particle colliders. 
Solution method: Simpson's rule integrator, basic algebraic computation. 
Additional comments: SPICE interfaces with SOFTSUSY and SUSYHIT to produce the low- 
energy sparticle spectrum. Flavor mixing for sleptons and sneutrinos is fully implemented; 
flavor mixing for squarks is not included. 



II. INTRODUCTION 

There are two flavor problems in particle physics. The standard model (SM) flavor 
problem is the unexplained pattern of SM fermion masses. In addition, there is the new 
physics flavor problem, that is, the difficulty of solving the gauge hierarchy problem, which 
inevitably requires the introduction of new particles at the weak scale, without generically 
violating well-known low-energy flavor constraints. Future colliders, including the Large 
Hadron Collider (LHC), may shed light on both the SM and the new physics flavor problems. 

Unfortunately, the potential of the LHC in this regard is not well explored, in part 
because simulations of new physics typically focus on the simplest possibilities in which 
flavor violation is absent. In this paper, we describe a new tool, SPICE: Simulation Package 
for Including Flavor in Collider Events, which is intended as a flrst step toward addressing 
this problem. 

In more detail, there are several motivations for developing collider event simulation tools 
for models with non-trivial flavor physics: 

• In many new physics frameworks, such as supersymmetry or extra dimensions, new 
states are predicted that are governed by the same flavor or horizontal symmetries as 
the SM fermions. Understanding their masses and mixings may therefore lead to real 
progress in identifying the theory of flavor. 

• Although the new physics may be minimally flavor violating, this is by no means 
required by current constraints; see, for example, Refs. j3, H, 0, 0, @, @, [13, 11 



As LHC data become available, it is appropriate to relax theoretical constraints and 
consider more general frameworks, especially if they yield unusual signals. 



2 



• The exploration of new physics is often said to consist of three phases: particle discov- 
ery, measurement of particle masses, and measurement of flavor and other mixings. In 
reality, these stages are unlikely to be completely distinct, and the exploration of flavor 
effects may provide new opportunities for particle discovery or proceed in parallel with 
precision mass measurements. 

• There are many promising probes of flavor violation at low energies. If a positive 
signal appears in one of them, for example, fi 67, there will be great interest 
in understanding whether new weak-scale physics can explain it. In this scenario, 
a program like SPICE will be helpful to understand what flavor- violating effects are 
observed or observable at colliders. 

• Finally, searches for new physics often implicitly assume the (near) absence of flavor 
violation. It is important to test their robustness in the presence of flavor violation. 

Although SPICE is flexible in several ways, in its typical form SPICE reads in parameters 
for one of the standard flavor-conserving supersymmetric models and lepton flavor- violating 
(LFV) masses specified by the user. SPICE uses the flavor-conserving model as a spine, 
then includes the corrections from LFV terms, determines the new superpartner masses and 
mass eigenstates, calculates branching ratios for all relevant flavor- conserving and flavor- 



violating decay modes, and outputs the results in HERWIG [12, ll3|, uM, uM SUSY Les 



Houches Accord (SLHA and SLHA2) jl6l. Il7l| format files, which may be used to generate 
collider events. SPICE uses SOFTSUSY |2| to evolve supersymmetric parameters through 



renormalization group equations and uses SUSYHIT [3[ to calculate flavor-conserving decay 
widths. 

The flavor- conserving model may be any of the standard ones, for example minimal 
supergravity (mSUGRA), minimal gauge- mediated supersymmetry breaking (mGMSB), 
or anomaly-mediated supersymmetry breaking supplemented by a universal scalar mass 
(mAMSB). We review these and their input parameters below. 

The LFV may be specified in a number of ways. At the "highest," most model-dependent 
level, a user may specify a horizontal global symmetry 18, 19, 20], the charges of each of 



the supermultiplets under this symmetry, and the Froggatt-Nielsen expansion parameter A. 
SPICE then populates each LFV mass matrix with entries that have the appropriate powers 
of A, multiplied by randomly-chosen 0{1) coefficients. At a slightly lower level, the user 
may specify each of the 0{1) coefficients by hand. And finally, at the lowest and most 
model-independent level, a user may specify every LFV term by hand. 

Supersymmetric models with LFV are, of course, just some of the many possibilities for 
non-trivial flavor physics at the LHC We are led to consider LFV in supersymmetric models 
for a variety of reasons. Supersymmetry is well-motivated by the gauge hierarchy problem 
and force unification, both with and without gravity, and there exist supersymmetric flavor 
models that can both explain the observed charged lepton and neutrino masses and mixings 
and simultaneously solve the new physics flavor problem. In addition, and perhaps most 
important for the coming data-driven era, supersymmetry is flexible enough to encompass 
many diverse experimental signatures, and lepton flavor violation, as opposed to quark flavor 
violation, generically has obvious implications for observable experimental signatures and 
event topologies. That said, it would certainly also be interesting to study hadronic flavor 
violation, particularly that involving the 3rd generation, and to generalize these results 
to non-supersymmetric frameworks. A recent study of the observability of squark flavor 



violation at the LHC was done in Ref. 21 
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The program and detailed installation instructions can be found on the web at the address 

http://hep.ps.uci.edu/~spice . 
In addition, App. |A] is a step-by-step guide to installing SPICE and gives helpful tips on 
getting started. 

In the following sections, we explain how SPICE works. In Sec. Illlt we discuss the the- 
oretical framework in more detail. In the presence of LFV, many new decay modes open 
up, but in most cases, the decay widths are simple generalizations of flavor- conserving ones. 
An exception is the case of "charge-flipping," three-body slepton decays to like-sign leptons 
i~ — > We have calculated these and discussed their interesting phenomenology 



elsewhere 22|] , but we summarize the main points in Sec. IIIII 

In Sec. llVt we show how the program works, explaining how to install the program and 
specifying the relevant file names and input and output formats. To aid first-time users, we 
also present a simple example input file and the resulting output. Our conclusions are given 
in Sec. El Further details regarding our conventions, Lagrangian terms, decay widths, and 
program details are given in a series of appendices. 



III. THEORETICAL FRAMEWORK 
A. Flavor- Conserving Inputs 

The user inputs to SPICE consist of a flavor-conserving SUSY spine supplemented by 
LFV terms. The flavor-conserving SUSY spine may be any model normally available in 
SOFTSUSY fll, 0] , though in principle any model could be input with appropriate alterations 
to the SOFTSUSY code. The base options and their input parameters are the following: 

• Minimal supergravity (mSUGRA) is specified by 

mSUGRA: mo,Mi/2, Ao,tan/3,mGUT,sign(/i) , (1) 

where mo, M1/2, and Aq are the universal scalar, gaugino, and tri-linear scalar coupling 
parameters at the grand unified theory (GUT) scale mcuT, tan/? = {Hq)/{Hq), and 
sign(yu) is the sign of the Higgsino mass parameter /i. 

• Minimal gauge-mediated supersymmetry breaking (mGMSB) is specified by 

mGMSB: N^, M^^ss, A, Cg^av, tan/?, sign(/i) , (2) 

where is the number of 5 + 5 multiplets in the messenger sector, Mmess is the 
messenger mass scale, A = Fs/{S) sets the mass scale for the SM superpartners, 
Cgrav = 1^0/ (F/^/Smpianck) is the gravltluo mass in units of its mass in the minimal 
case in which there is only one SUSY-breaking F term, and tan/? and sign(/i) are as 
above. 

• Minimal anomaly-mediated supersymmetry breaking (mAMSB) is specified by 

mAMSB: mo, ^3/2, tan/5, mcuT, sign(/i) , (3) 

where mo is the universal scalar mass, motivated by the tachyonic slepton problem 
and added to all scalars, including the Higgs scalars, m3/2 is the gravitino mass, which 
sets the scale for all SM superpartners, and tan/5, mcuT, and sign(/i) are as above. 
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B. Flavor- Violating Inputs: Model-Independent Approach 

The LFV terms are specified by a parameter x and three 3x3 matrices itle, Xl, and X^. 
All of these matrices are assumed real, and X^ and Xr are necessarily symmetric. In the 
most model-independent formulation accommodated by SPICE, LFV is therefore completely 
specified by 22 numbers: x, 9 for rriE, 6 for X^, and 6 for Xr. 

The matrix m^; is the SM charged lepton mass matrix. X^ and X^ determine the slepton 
mass matrices through 

M? = mll + xm^XL (4) 
= mll + mEml; + xm^XL (5) 
— "^^1 + 'ni^E'^E + xfh^Xji , (6) 

where m| and m| are the flavor-conserving contributions to the left- and right-handed 
sleptons, and characterizes the size of flavor-conserving contributions to the slepton 
masses. For concreteness we take to be equal to the average of the diagonal elements of 
m| ^. The parameter x specifies the size of the LFV effects relative to the flavor-conserving 
parameters. 

Several simplifying assumptions are encoded in this formulation of the LFV effects: 

1. In assuming mg, X^, and Xr are real, we do not include CP violation from these 
matrices. 

2. The flavor- violating contributions to left- and right-handed sleptons are of the same 
order, as set by xin^. This is what one would expect of gravitational contributions, 
which are chirality-blind. 

3. In general, there should also be LFV A-terms. We assume these are negligible at the 
mass scale where these SUSY-breaking masses are generated. Note, however, that 
^-terms are generated through RG evolution, and such effects are included. 

4. In general, there should also be a non-trivial LFV neutrino mass matrix. However, 
since we are primarily interested in colliders, where neutrino flavor is unobservable, we 
simply assume that the neutrino masses vanish, and so the neutrino mass and gauge 
eigenstates are identical. Note that sneutrinos may have observable mixings, and these 
are included. 



C. Flavor- Violating Inputs: Flavor Symmetry Approach 



An attractive possibility is that the LFV terms are determined by horizontal symmetries. 
In this approach, the 3 LFV matrices are 



ruE = me. 



I ciX^^ C2A"2 csA"^ 
C4A"* cgA"^ ceA"^ 

^C7A"7 CgA"* CgA'*^ 



(7) 



^ For a flavor-conserving spine where the slepton mass-squared matrices are not proportional to the identity, 
w? is thus equal to the average left slepton mass-squared. 



5 



Xt 



X 



R 



cioA^io ciiA"ii ciaA'^i^ 

ciiA"" C13A"" C14A"" 

Cl2A"i2 C14A"" Ci5A"i5 

/cigA"" citA"^^ ci8A"i« 

cirA'^i^ C19A"" C2oA"2o 

VCi8A"i« C2oA"2« C2iA"2i 



(8) 
(9) 



where is the lepton mass scale, the Cj are 0{1) coefficients, A is the Froggatt-Nielsen 
expansion parameter, and the exponents are determined by supermuhiplet charges. 

SPICE accommodates this possibihty by providing an alternative way to specify the LFV 
parameters in the case of [/(l)" flavor models, with breaking parameters of the same size. 
For example, for U{1) x f/(l), one sets A^'charges = 2, and specifies the two U(l) charges for 
the six multiplets Li, L2, L3, Ei, E2, and E^. SPICE then determines the correct exponent 
for each entry of the three matrices. SPICE may also determine the coefficients q randomly 
or these may be set by the user. 



D. Spectrum and Decay Calculations 

Once the flavor-conserving and flavor- violating parameters are input, the model is com- 
pletely specified. Formally, one should then evolve the mass parameters from the scale 
they are generated, adding in new contributions at the appropriate scale, until one reaches 
the weak scale. In SPICE, we instead add the LFV contributions to the flavor- conserving 
contributions at the scale at which the flavor-conserving contributions are generated, and 
then RG evolve the result to the weak scale. In the case that the flavor- conserving and 
LFV contributions are generated at the same scale, as in the case of mSUGRA or mAMSB 
spines with LFV generated at mQUT, this is the correct prescription. If the generation scales 
differ, for example, as in the case of a mCMSB spine with gravity-mediated LFV effects, 
our prescription is not formally correct. Even in this case, however, we expect that the 
dominant difference resulting from this approximation may be absorbed into the LFV 0{1) 
parameters, which are not completely specified anyway. 

Given the slepton masses and mixings determined above, we then calculate all branching 
ratios involving these LFV effects. Most of the flavor-violating decay widths are, at least to 
tree level, fairly simple generalizations of the flavor-conserving widths. New effects appear 
in the 3-body decays of sleptons both to like-sign leptons, — *• i'j'ik^rny opposite- 



sign leptons, i~ iji^i^. Previous work on these decays [23|, |2J| has assumed flavor 
conservation. In the presence of LFV, however, the lepton pair may be identical particles 
or charge conjugates of the same generation, leading to additional interference terms. We 
have calculated these decay widths in the presence of LFV and arbitrary left-right slepton 



cay 

mixing in another work 22| and included these in SPICE. 

Currently SPICE discards decays to gravitinos. Such decays are typically either kinemat- 
ically inaccessible or too slow to be relevant at colliders. Even in the cases where the model 
used for the spine predicts fast decays to gravitinos, as in the case of mGMSB with a low 
supersymmetry-breaking scale, the existence of significant LFV contributions implies that 
in the full LFV theory, supersymmetry is broken at a high scale, and decays to gravitinos 
are negligible. 
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1. User-supplied SPICEinit file is read 
by spice.x. 

2. spice.x uses SOI=TSUSY to generate 
the mass spectrum and calculates 
lepton flavor-violating decay widths. 

3. SUSYHIT uses output from spice.x 
and calculates flavor-conserving 
decay widths. 

4. FileCombine.x combines the output 
from spice.x and SUSYHrT to create 
HenwigFinal.out, SLHAFinal.out, and 
SLHAZFinal.out, which can be used 

to generate collider events. 



SPICEinit 



spice.x (SOFTSUSY) 



softsusySLHA.out 

renamed to 
slhaspectrum.in 









SUSYHIT 






susyhit_slha.out 



spicydecays.out 
softsu sySLH A2 .out 



FileCombine.x 



HerwigFlnal.out 
SLHAFinal.out 
SLHAZFinal.out 



FIG. 1: Schematic flow chart for SPICE, showing the progression from the SPICEinit input file 
to the final output files HerwigFinal . out, SLHAFinal.out, and SLHA2Final . out. The explicit 
interfacing with SOFTSUSY Q, [J and SUSYHIT jj^] is also shown. 



IV. PROGRAM 
A. Procedure 

The program is designed to take in one input file which specifies all of the input param- 
eters. This file, by default, is called SPICEinit and is located in the spice/ directory. The 
default input file is identical to the example input file presented below. These input param- 
eters are used by SOFTSUSY 0] to generate a particle mass spectrum, which is passed to 
SUSYHIT 1^ to create a SUSY Les Houches Accord (SLHA) file with the traditional flavor- 
conserving decay table. Decay widths of the flavor-generalized decay channels involving 
sleptons and sneutrinos are calculated separately within the SOFTSUSY executable, which 
generates an intermediate file containing the flavor-generalized mass spectrum and relevant 
decays. The FileCombine subprogram merges the SUSYHIT output with this intermediate 
file to generate the file HerwigFinal . out, which is an input file for the Herwig event gen- 
erator [ii, EEJli, [111 and both SLHAFinal . out and SLHA2Final . out, which are SLHA ^ 



and SLHA2 [17[ formatted files appropriate for input to Monte Carlo event generators. The 
program flow is diagrammed schematically in Fig. [TJ 

The input file for the program contains two sets of input parameters — the flavor- 
conserving SUSY spine, and the flavor-violating effects. The first set specifies the SUSY 
breaking scenario parameters according to the SOFTSUSY [l], 0| specifications, which are 
briefly reviewed below. The second set describe high energy flavor mixing boundary condi- 
tions. The default input file is presented below as an example^. 



This example file uses the U(l) horizontal charges of Model B in Ref. 
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Example Input File 



gmsb 4 2.0e6 5.0e4 1.0 10 1 
X 0.1 

lambda 0.2 
nCharges 2 
LI 2 
L2 2 
L3 2 
El 2 1 
E2 2 -1 
E3 -1 

Lep -0.13854 2.24310 -0.4399 1.84877 1.53904 -0.56145 0.89942 0.53290 
1.45314 

XL 0.98914 -1.2001 1.40588 -1.2001 -0.470473 2.67457 1.40588 2.67457 
-0.22933 

XR -0.5889 2.8415 -0.2321 2.8415 -0.43167 1.3267 -0.2321 1.3267 0.55438 

• Flavor conserving parameters: 

We incorporate tlie tliree SUSY breaking scenarios implemented by SOFTSUSY as 
options for the flavor-conserving SUSY spine. These are mSUGRA, mGMSB, and 
mAMSB, as described in Sec. Illll The flrst line of the script flle specifles which spine 
is used with an identifler of "sugra" , "gmsb" , or "amsb" , and then parameters needed 
for that given scenario. The ordering is consistent with that used in Sec. Illll and the 
SOFTSUSY manual [l|, 0]. Note that the text labels for the SUSY breaking scenario is 
case sensitive. These are summarized as follows: 

sugra <mO> <ml2> <aO> <tanb> <mgut> <sgnMu> 

gmsb <n5> <mMess> <lambda> <cgrav> <tanb> <sgnMu> 

amsb <mO> <m32> <tanb> <mgut> <sgnMu> 

• Flavor mixing parameters 

The remainder of the input flle is used to describe the slepton flavor mixing in the 
model in the context of U(l) horizontal charges as described in Ref. jl] and reviewed in 
Sec. Illll The flrst two parameters are x and A, and the next integer nCharges specifles 
the number of U(l) charges for the leptons. Note that text labels are required as 
placeholders before each variable, although the actual label text is arbitrary for all 
entries except the spine identifler and the "Lep" identifler which may be changed to 
"random" to generate random coefficients (see below). After the number of charges is 
passed, the charge assignments for each generation of left-handed sleptons and right- 
handed sleptons must be given. Lastly, the user can specify the 0{1) coefficients to be 
used in the lepton, left-handed slepton, and right-handed slepton mass matrices. The 
blocks "Lep", "SLepXL", and "SLepXR" define the nine 0{1) coefficients for each 
3x3 mass matrix, according to the pattern 



Lep Ci C2 C3 C4 C5 Ce C7 Cs Cg 



C4A"* csA"^ C6A"6 

\C7A"^ CgA"* CgA"^ 



) 



(10) 
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"SLepXL" and "SLepXR" are symmetric matrices and must be defined with symmet- 
ric coefficients, while the "Lep" matrix has arbitrary coefficients. Alternatively, the 
randomly generated coefficients may be generated by replacing the last three lines with 

random <rseed> <sigma> 
where <rseed> is an integer, used as a random seed to generate 0{1) numbers of 
the form ±exp(a), where a has a gaussian distribution with a mean of and a a of 
<sigma>. These coefficients have a 50% probability of being negative. 

Besides the default example file, there are several example input files provided in the sub- 
directory spice/Examples/ that demonstrate various charge assignments following the mod- 
els A-D prescription in Ref. To use these, rename them to the default input file name 
SPICEinit in the spice/ directory, or modify the CalcDecays command script appropri- 
ately. 

B. SOFTSUSY and SUSYHIT Details 

SPICE relies upon SOFTSUSY jH, [2|, using SOFTSUSY for renormalization with flavor-general 
boundary conditions. The lepton flavor mixing is rotated into slepton and sneutrino mass 
matrices at the high scale before renormalization is performed; the rotation of the neutrino 
masses is neglected along with the neutrino masses themselves. In general, diagonalizing 
the lepton mass matrix requires rotating the trilinear slepton couplings as well, but this is 
unnecessary in our case as we assume negligible trilinear contributions at the high scale. 
Using the flavor general boundary conditions in spice. cpp, SOFTSUSY iteratively solves the 
renormalization group equations (RGEs) and outputs a SUSY mass hierarchy in SLHA and 
SLHA2 format to sof tsusySLHA . out and sof tsusySLHA2 . out, respectively; the SLHA ffie 
is passed to SUSYHIT as slhaspectrum. in. We do not force gauge coupling unification, the 
GUT scale is set at mcuT = 2x10^^ GeV initially, we use 2-loop Higgs formulas for physical 
masses, and we generally follow SOFTSUSY program defaults and conventions as much as 
possible. For more details on the workings of SOFTSUSY, refer to the SOFTSUSY manual 0] . 

SUSYHIT takes slhaspectrum. in and calculates non-flavor violating decays of the entire 
SUSY spectrum, according to its subroutines SDECAY and HDECAY. Details can be found in 
the SUSYHIT reference manul (sf. The output file appends the calculated decay table to the 
input ffie slhaspectrum. in, and this SUSYHIT output ffie is passed to FileCombine. 

C. Decay Calculation Files 

The calculations of new flavor violating decays is performed within the SPICE 
executable with the code contained in DecayFormat . cpp, DecayCalc . cpp and 
ThreeBodylntegrals . cpp and their respective .h header ffies. The content of these ffies is: 

• The DecayFormat class: 

— Output of particle mass and lifetime data 

— Output of a standard format of two- and three-body decays 

— Output of SUSY parameters in HERWIG format 
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— A public method which chains these methods together to produce a pseudo- 
HERWIG file^ 

• The DecayCalc class: 

— Retrieval of the slepton/sneutrino mass matrices from FlavourMssmSof tsusy 

— Calculations for all decays involving sleptons and sneutrinos 

— Functions for Lagrangian coefficients and three-body decay integral multiplicative 
factors 

— Loop over all possible sleptons/leptons for the given decays and output via 
DecayFormat class 

• The ThreeBodylntegrals class: 

— Functions for the integrands of the three body decay integrals 

— A Simpson's rule integration method 

— Methods to calculate the normal and interference term integrals given a vector 
containing multiplicative factors 

The decay calculations (20 in total) are numbered according to Table II in App. [Dl 
These replace all the standard decays calculated by SUSYHIT that involve sleptons/sneutrinos 
as either decaying particle or product. The routine is called from spice, cpp using the 
OutputHerwigFile method of a DecayFormat object, which in turn creates a DecayCalc 
object and proceeds to call the other global output functions to generate a file in pseudo- 
HERWIG format containing the slepton/sneutrino decays. This method also generates the 
mass spectrum and model parameters for inclusion in the final output file appropriate for 
Monte Carlo event generators. There are also several decay modes listed in Table III in 
App. [D] that are irrelevant or not viable in the models we have studied but may be viable 
in other scenarios; instructions on how to insert new decay modes are presented in App. O 

D. Merging Flavor Violation and SUSYHIT Files 

The final component of SPICE is the FileCombine subroutine, which takes the fiavor- 
violating slepton/sneutrino decays file and merges it with the other decays calculated 
by SUSYHIT. FileCombine takes the particle mass info and model parameters from the 
softsusySLHA.out file, along with the slepton-sneutrino decay table. It then reads through 
the SUSYHIT decay file and stores all decays that do not involve sleptons or sneutrinos. 



^ This pseudo-HERWIG file differs from a true HERWIG file since it lacks SUSY decays beyond those involving 
sleptons and sneutrinos, it does not have the appropriate whitespace formatting, and it uses Particle Data 
Group (PDG) particle identification codes instead of HERWIG particle codes. 
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FIG. 2: Low-energy spectrum and decay modes for the flavor-conserving GMSB spine specified 
in the SPICEinit example input file. Flavor-violating effects were not included, and decays to 
gravitino are ignored. For the slepton shadings (colors), stau flavor is medium gray (red), smuon 
flavor is light gray (yellow), and selectron flavor is dark gray (blue). SUSY decays proceed along 
arrows and are labeled by the associated SM daughter particles. The relevant masses are x? = 



266.7 GeV, 



200.1 GeV, ^2 = 200.1 GeV, and 



193.6 GeV. 



E. Output 



The program outputs the stored data as a standard HERWIG input file 1^, 1^, 14, 15] and 
both SLHA- and SLHA2-formatted files for input to other Monte Carlo event simulation 
packages. The first portion of the HERWIG file contains particle codes, masses, and lifetimes; 
the bulk contains branching fractions for particle decays, and the final section contains 
model parameters^. It contains all feasible (and large) flavor- violating slepton and sneutrino 
decays. 



For the example input file in Subsection IIV Al we provide some comments about the expected 
output and how our example model differs from traditional flavor-conserving models. Using 



The SLHA files are structured similarly; for details, see Refs. 
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FIG. 3: Low-energy spectrum and decay modes for the full flavor-violating model specified in 
the SPICEinit example input file. Decays to gravitino are ignored. Flavor-violating effects are 
included and serve to modify slepton masses, mix the gauge eigenstates, and introduce many new 
decay modes. For the slepton shadings (colors), stau flavor is medium gray (red), smuon flavor is 
light gray (yellow), and selectron flavor is dark gray (blue). SUSY decays proceed along arrows 
and are labeled by the associated SM daughter particles. The relevant masses are Xi = 266.7 GeV, 
q = 210.0 GeV, q = 188.7 GeV, and i{ = 179.8 GeV, and the branching ratios are listed in 
Table H 

only the flavor- conserving spine of the example model, for instance, we derive the traditional 
low-energy spectrum and straightforward decay channels depicted in Fig. [2J When SPICE 
calculates the full flavor-violating example model, however, the decays modes acquire new 
complexity and richness: these modes are diagrammed in Fig. [3] and listed in Table [B We 
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7.7% 


6.5% 








4 more 
0.3% 




210.0 




ife-T- 
44.3% 


23.9% 


9.5% 


8.4% 


5.2% 


4.8% 


ete-T- 

1.9% 


1.1% 


50 more 
0.9% 


^2 


188.7 




27.8% 


21.6% 


ife'e- 
20.0% 


19.7% 


4.7% 


4.7% 


1.0% 




22 more 
0.5% 



TABLE I: List of masses (in GeV) and branching ratios for decaying particles Xi; ^3 ; ^^nd £2 using 
the SPICEinit example input file. Only branching ratios larger than 1.0% are individually listed, 
and decays to gravitino are ignored. This table refers to the spectrum depicted in Fig. [3l 
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note that the flavor mixing parameters for the Xr matrix in the example input flle indicate 
that the two hghtest sleptons are predoniinantly selectron and smuon, which is conflrmed 
by the branching ratios for Xi ~^ ifi^,ifi^. Thus, while traditional expectations favor 
the lightest slepton to be predominantly stau, given the large Yukawa coupling of the third 
generation, the lightest sleptons can instead generically be combinations of selectron and 
smuon in flavor-violating scenarios. More importantly, the wealth of new information in 
flavor-violating two- and three-body modes indicates it may be possible to measure flavor 
couplings from early collider data, if positive signals for SUSY are discovered. 



V. CONCLUSION 

We have presented the program SPICE: Simulation Package for Including Flavor in Col- 
lider Events intended for use in studies of lepton flavor violation (LFV). SPICE simulates 
hybrid models where flavor-violating gravitational effects are added to a flavor-conserving 
SUSY spine. The program contains preset options for generating lepton flavor violation 
based on U(l) horizontal symmetries, but it also allows for input of more general flavor- 
violating effects. It interfaces with SOFTSUSY, which generates the low-energy mass hierarchy, 
and SUSYHIT, which calculates decays unconnected to LFV. SPICE automatically calculates 
all kinematically allowed decays involving sleptons and sneutrinos with LFV, and also cal- 
culates three-body slepton decays, which are important in models with a long-lived slepton. 
Given an input model with LFV, SPICE outputs the mass hierarchy and decay table with 
LFV effects, which can be used to generate collider events using Monte Carlo simulation 
packages. SPICE is freely available under the Gnu Public License, and can be downloaded 
at http://hep.ps.uci.edu/~spice . 

There are several possible extensions to SPICE. Foremost among these is the inclusion 
of quark sector flavor violation, which, while not expected to impact collider signals as 
obviously as some possible LFV signals, also provides interesting effects. Another valuable 
extension would be the inclusion of low-energy effects, including neutrino masses and mixings 
and comparison to low-energy constraints such as /i — 67. Another possibility is the flavor 



generalization and inclusion of three-body decays involving sneutrinos as detailed in Ref. [24 
In the case of a long-lived slepton, there is also the possibility of a decay of the form ii £^7 
generated through LFV loops, which could be considered. 
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APPENDICES 



This section of the paper provides all details relevant to running SPICE. Appendix Rl 
outlines the installation procedure and provides general tips on getting started, as well as 
basic troubleshooting solutions. Detailed command-line instructions can be found in the 
README file on the website. Appendix [B] covers the notation and conventions used in SPICE. 
Appendix O defines the relevant portions of the Lagrangian for fiavor-general sleptons and 
defines Lagrangian coefficients to simplify the forms of the decay modes. Appendix [D] enu- 
merates the flavor- generahzed decays used and calculated by SPICE, as well as the decays 
SPICE ignores. Appendix [E] gives formulas for flavor-generalized two-body decays. Ap- 
pendix [F] discusses the three-body decays; the full formulas are presented in Ref . |22| . Ap- 
pendix [Gl flnally, discusses the structure of the program and details about the integration 
between SPICE, SOFTSUSY, and SUSYHIT. 

APPENDIX A: INSTALLATION INSTRUCTIONS AND TROUBLESHOOTING 

1. Installing and Running SPICE 

SPICE can be downloaded from 'http : / /hep . ps . uci . edu/ ~spi ce I . SPICE 
requires both SOFTSUSY (http://projects.hepforge.org/softsusy) and SUSYHIT 
(http://lappweb.in2p3.fr/~muehlleitner/SUSY-HITj) to be installed in order to func- 
tion properly. After obtaining and unzipping all three packages, use the SPICE Makefile 
in the spice/ directory to build the executables spice. x and FileCombine .x, and also 
run the default SUSYHIT Makefile to compile SUSYHIT^. The command script CalcDecays, 
located in the spice/ directory, handles program execution and movement of input and out- 
put flies: by default, it uses the flle SPICEinit for input parameters. The output of SPICE 
is a Herwig input flle HerwigFinal . out, a SUSY LesHouches Accord flle SLHAFinal . out, 
and a SLHA2 flle SLHA2Final . out, all located in spice/FileCombine/. When running the 
default SPICEinit flle, the program output should match the example output flies in the 
spice/FileCombine/ directory. Additional input flies, which demonstrate the versatility of 
SPICE, can be found in spice/Examples/. 

2. Troubleshooting 

By default, SPICE assumes a particular directory structure, with the spice/, softsusy/, 
and susyhit/ directories all located in the same parent directory, i.e., Desktop/spice/, 
Desktop/sof tsusy/, and Desktop/susyhit/. Any other directory structure or different 
directory names can be accomodated by making the appropriate changes to the directory 
path names in the SPICE Makefile and command script CalcDecays. 

SPICE was tested with the standard SOFTSUSY 3.0.2. and SUSYHIT 1.3 builds, so a cus- 
tomized build of either program's core flies may cause errors. In addition, SPICE was com- 
piled using the g++-3 . 4 compiler, so use of a different compiler may require changes to 



^ Please note that spice. x contains the relevant portions of SOFTSUSY, so a separate SOFTSUSY build is not 
required. 



14 



the SPICE Makefile. In particular, for the gcc-4.2 compiler, the user should remove the 
-f f 90 option in the Makefile and replace the -lg2c inline tag with -Igf ortran. Further 
installation details are provided in the README file included with SPICE. 



APPENDIX B: CONVENTIONS 

In this section we use the spacetime metric (H ) and fermionic propagators of the 

form — rnj / (p^ — m^). Our projection operators are Pl,r — (1 ^75)72. The constants g 

and g' are the standard weak scale SU(2) and U(l) gauge couplings; we break from SOFTSUSY 
convention and use the typical standard model value for g' rather than the gauge unification 
value. 

Throughout this appendix, greek indices refer to gauge eigenstates while roman indices 
refer to mass eigenstates for clarity. This convention is dropped in later sections to avoid 
confusion with Lorentz indices. 

1. Higgs Bosons 

Our convention for the Higgs doublets is 




(Bl) 



and we set their vacuum expectation values 



(VEVs) to be 




(B2) 



We set the SM Higgs VEV v to be v'^ = vl + v\ 
spondingly, the SM gauge boson masses are 



and also we define tan/5 = v^jv^. Corre- 





(B3) 



and the SM lepton masses are 




(B4) 



After electroweak symmetry breaking, the neutral Higgs doublet is given by 




^0 



) 



(B5) 




(B6) 



where 




and 




(B7) 
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2. Neutralinos 

The neutralino gauge eigenstates are defined by = (^—iB,—iW,^p'^,^^"^ , with the 
ordering of the up/down higgsinos chosen to conform with SOFTSUSY. All the above fields 
are two-component Weyl spinors. 

The neutralino mass term in the Lagrangian is given by 

jC = _iv;°tM^o,^^2 , (B8) 

so diagonalizing according to Oj^M^o^^Opj = m^oSij gives neutralino mass eigenstates of 

Xi = Oj^t/j'^. The corresponding Dirac spinors are x° = (x^ Here we follow SOFTSUSY 
convention for the neutralino mass matrix by keeping the matrix real and allowing the 
neutralino masses to be negative. 



3. Charginos 

Similarly, the chargino gauge eigenstates are defined by V^"*" = (^—iW~^,tp:^^ and tp~ — 

(^—iW~,'tp'^^ . Again, these are two-component Weyl spinors. 
The Lagrangian mass term is given by 

so diagonalizing with the rotation angles 9l and 9r gives 

/ cos^L sin^iA^ f cosOr -smeR\ ^ I \ 

V-sin^i cos^l; '^^l^sin^R cos^r/ ^ j ' ^ ' 

The chargino mass eigenstates thus become 

ixt\^i cos Qr sin Qr \ ( -tW+ \ .g^^. 
\xt) K-sinOR cosOr) \ i^t ) 

and 

(Xi\^( cos Ol sin ~iW~ \ ^g^2) 

V X2 / V - sin Ol cos Ol) \ / 

The corresponding Dirac spinors are Xi — {xt iXi^ ■ 



4. Charged Sleptons 



The mass term for charged sleptons is 

Here la = {^l, P'L,t^l,^r, P'RjT^r)^ is a column vector containing the left-handed sleptons in 
the first three components and the right-handed sleptons in the last three components. 
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The 6x6 mass- squared matrix for sleptons is 



The mass matrix is diagonahzed according to {Mf^ ^ = uj^i (^|) .. ^^/3''^ so the slepton 

mass eigenstates are defined by ii = U-^^ia and = U^iii-, with the mass eigenstates 
ordered in terms of increasing mass. 

Note that the first index of the matrix t/ is a gauge index and the second is a mass index; 
almost all interactions treat the left- and right- handed sleptons differently, so typically the 
first index will be summed from either 1 to 3 or 4 to 6, while the second will be summed 
over all six lepton mass eigenstates. Note also that in our convention the rotation matrix is 



the hermitian conjugate of the rotation matrix presented in Ref. jl7| . 



5. Sneutrinos 



The sneutrino mass term is 

£ = z>:(m^)^^z>^, (B15) 

with = {i'e.v^^i'rf . Diagonaliziug gives (M|)„^ = U^ai\M?)-^U-p^\ and thus i>i = 

l&^i'a and Ua = Uj^j^Ui. Again our rotation matrix is the hermitian conjugate of the one in 
17 . 



6. Leptons 

As previously mentioned, we work in the basis of charged leptons in which the mass and 
fiavor eigenstates are in correspondence. This is possible because there is a freedom of a 
single superfield rotation, and we choose to use this rotation to rotate away any mixing in the 
charged lepton sector (this was done for both notational simplicity and convenient coding). 
Also, as mentioned previously, we have neglected neutrino mixing, as neutrino fiavor is not 
observable at colliders. 



APPENDIX C: LAGRANGIANS 

Indices in Lagrangian coefficients are ordered by "gaugino, sfermion, fermion" if there is 
a gaugino present. Our Lagrangian assumes no neutrino masses and thus that the charged 
leptons can freely be rotated into mass eigenstates. Thus all the mixing terms will be 
parameterized by only the slepton and sneutrino mixing matrices (and the standard neu- 
tralino/chargino mixing matrices). 

For the mixing matrices, the left and right bases are different, with one being the gauge 
basis and the other the mass basis. We have chosen the first index to correspond to a gauge 

eigenstate and the second to correspond to a mass eigenstate. Thus, in a term like U^2-> ^ 
the gauge index and 6 is a slepton mass index. Note that in summations over these indices, 
the gauge summations are always taken from 1 to 3 and the mass summations are taken 
from 1 to 6 for sleptons and 1 to 3 for sneutrinos. 
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In this section, we only consider phenomenologically viable decays to lowest order; thus 
we neglect four scalar couplings and couplings for modes that are unlikely to be kinematically 
allowed. 



1. Neutralino Terms 

The neutralino interaction terms are 

>CNeutralino = [aabci^^xiPLl^c + (/^S^Pl + /SSi^fi) ic] + h.C. , (CI) 



where 



C^abc = {gOla - g'Ola) t^if (C2) 

/3« = i= {gOl^ + g'Ol^) U$; - y\^^OlJJ% (C3) 
& - -V20,,aM% - y^'0,,M>. (C4) 

For these terms, the matrix Oij is the neutralino mixing matrix and y^^ is the c component 
of the diagonal lepton Yukawa matrix; the gauge components correspond to 1 - bino, 2 - 
wino, 3 - down Higgsino^. 

2. Chargino Terms 

The chargino interaction terms are 

>Cchaxgino = {labMaPLl^c + ^t¥a {^^^Pl + ^Si^s) 4} + h.C, (C5) 



where 



Tabc 



-g cos Ol Uf}* + sin Ol C/i+a* 6 « = 1 /^gN 
gsine,uf> + y^)cose,U% a = 2 



(1) ^ \ -g cos 9r u''^^* a = 1 
[ g sm Or U^^ a = 2 



(C7) 



a(2) _ J #sm^Lf/y a = l , . 

^abc ~ \ (£) n TTi'^)* n ' ^ ^ 

[ y^^> cos a*!, U'^^^' a = 2 

Here and 9r are the rotation angles associated with rotating the negatively and positively 
charged (respectively) gauginos and higginos into the mass eigenstates. 



The up Higgsino corresponds to o = 4, but this is never needed. 



18 



3. Higgs Terms 



The neutral Higgs interaction terms are 

/:Neutrai Higgs = cy^^ K^H' + ^i^C^/^" + alJc^i/" + la^^tM' + ia^^l%G\ (C9) 

There are five different a coefficients since the two slepton couphngs to and are 
different: 



a 



(1) 

ah 



gmw 



a 



(2) 
ah 



2 cos^ 9w 
( gmw 
V 2 



cos(6'h + j3)5ah 

1 - tan^ Ow) sin(6'H + /?) 



(CIO) 



gfm^^sin Oh 
mw cos (3 



+ gmw tan^ 6'^/ sin(6'iy + /?) — I U!:Ar.U. 



cr 



(3) 
ab 



2mw cos /? 
I — ^ — (1 — tan 



mw cos P 

(/.cos^H + sin^^) (f/(2*f/ij3„ + t/iit^S 

(yfmf^ COS 6h 



(Cll) 



+ I (7mpt/ tan^ 6^1^ cos{6h + /?) 
4 



mvi^ cos P 
gmj cosOh' 



'-^ c+3,a'-^ c+3,1 



mw COS P 



2mw cos P 



c,h 



a 



(4) 
ah 



2mw 

gm^ 

2mw 



(/itan/3 + Af)U/(2*f/i 



(C12) 
(C13) 
(C14) 



Here tan^vv = g'/g as usual, tan/5 = Vu/vd, and the angle 9h is the rotation angle for the 
neutral Higgs scalars Hu and Ha into the real neutral Higgs mass eigenstates. The light and 
heavy real Higgs bosons are denoted and H^, and the pseudoscalar Higgs is denoted A^. 
The parameter /i is the Higgsino mass parameter, and Ai- is the left-right mixing term for 
ii in the diagonal lepton basis. 

It is useful to note that several of the neutral Higgs terms are proportional to U^^^*Uc^bi 

Uc+3,aUc% hi or uf}.l^Jjfl. The first two are "block fiavor diagonal" in left and right sleptons 
and the third has explicit left-right mixing, so fiavor mixing decays of these modes should 
be significantly suppressed relative to chargino/neutralino modes by the size of left-right 
mixing, as discussed in Ref. [2^ . 

Using the conventions above, the charged Higgs interaction terms are 



where 



Charged Higgs 



Pah 



+ h.c. 



(C15) 



mw 



mj tan/3\ (i\ 
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(2) 

Pab = 9 



( nT-w 



{l, + AeJanf3)Ui^*U% 



V2 



cos 2/3 - 



(C16) 



-^(/.tan/3-.l,jC/ifC/i53,, 



c,a c,b 
(1) 



(C17) 



Note there is no additional rotation angle since /3 serves as the rotation angle for the charged 
Higgs. 



4. Gauge Bosons 

The Lagrangian interaction terms for the gauge bosons are 



^Gauge Boson 



Ki^ iKd,h - ibd.K) + h.c.; 
+iCi? {ildj, - hdJi) + icif {W, - hdJi) , (C18) 



where 



Safe 

Safe 
.(3) 



3_TT{i>)*Tjii) 

^^^c,a ^c,fe 



9 



2 cos 9w 



(C19) 

(C20) 
(C21) 



Here, e is the electromagnetic charge. As above for neutral Higgs bosons, the flavor mixing 
modes involving Z bosons are heavily suppressed by the left-right mixing. The photon vertex 
remains flavor diagonal even in the flavor mixing case, as expected. 



20 



APPENDIX D: TABLE OF DECAY MODES 



Decay Number 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



Decay Mode 



Xi 
Xi 



70 



Vi - 
Vi - 

H- 
l7 - 



xjf-k 

' ij^k^m 



Corresponding Program Method 



CalcChiOSlepLepBar 

CalcChiOSnuNuBar 

CalcChiMinusSlepNuBar 

CalcChiPlusSnuLepBar 

CalcSlepChiOLep 

CalcSlepChiMinusNu 

CalcSlepSlepHiggs 

CalcSlepSlepZ 

CalcSlepSnuWMinus 

CalcSnuChiONu 

CalcSnuChiPlusLep 

CalcSnuSlepWPlus 

CalcHSlepSlepBar 

CalcHSnuSnuBar 

CalcASlepSlepBar 

CalcHMinusSlepSnuBar 

CalcSlepSlepLepLepBar 

CalcSlepSlepBarLepLep 

CalcSlepSlepNuNuBar 

CalcSlepSlepQQBar 



Table II: List of included decay modes with corresponding program methods in order of 

appearance. 



Neglected Decay Modes 



Justification 



^7 

- 
- 

i>i - 
h - 



P-H+ 



Any with emitted 7 
Scalar three-body decays 



assumed heavier than £ 
assumed heavier than i 
H~ assumed heavier than I 
Z assumed heavier than z> splitting 
hP assumed heavier than u splitting 
assumed heavier than u 
assumed heavier than v 
multiple £'s assumed heavier than 
multiple i/'s assumed heavier than hP 
QED interactions still flavor neutral 
Suppression by couplings and phase space factors 



Kinematics 
Kinematics 
Kinematics 
Kinematics 
Kinematics 
Kinematics 
Kinematics 
Kinematics 
Kinematics 



Table III: List of neglected decay modes, and reason for omission. 
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APPENDIX E: TWO BODY DECAYS 



This section contains the matrix elements and decay widths for all two-body decays. 
Charge conjugate modes are not included; we work in a non-CP violating lepton sector, 
so charge conjugate modes in all cases have the same decay widths as the original modes. 
Again, neutrinos are taken to be massless. 

Our index notation is that the index i refers to the decaying particle and the indices j and 
k to the decay products. Further, if the decaying particle is a slepton, sneutrino, chargino, 
or neutrahno then the index j refers to the resultant SUSY particle while k refers to the 
resultant Standard Model particle; this convention does not apply to Higgs decays, which 
have two SUSY daughters. 

The general formula for the decay width of a two-body decay in the center of mass frame 
of the decaying i particle with daughters j and k is 

I Alp 

r2 body decay = ^6^^^^^ ' ^^^^ 

where X{x, y, z) = + if + — 2{xy + xz + yz), and the matrix element squared has an 
implicit average/sum over fermion spins if the decay involves fermions. For our purposes, 
there are three general cases of matrix element forms: those containing two fermions (lepton 
and neutralino/chargino), those containing gauge bosons, and those containing only sleptons 
and Higgs scalars. The scalar case is trivial, since the matrix element is equal to the 
Lagrangian term with no momentum dependence, and the decay width is simply the square 
of the Lagragian term multiplied by a phase space factor. For fermionic modes with fermions 
/i and /2 and a scalar s, the matrix element takes the form of 

Atfe.n.ionic 2-body ^ ^ ( ^.f'^ ] {APn + BP,) ( ) , (E2) 

which gives a decay width proportional to ±{\A\'^+\B\'^){m'j_^+mj:^—'m'l)±4:mf-^rnf2 Re (AB*) 
times the standard phase space factor (with the possibility of m/j = if a neutrino is a decay 
product) . The sign on the second term will be positive if both spinors are uorv, and negative 
if they are different spinors. 

Finally, gauge boson matrix elements with scalars Si and S2 and a gauge boson b are pro- 
portional to {psi +Ps2Y^iJ-i so the decay width will be proportional to A2 (mg^, m^^, j /m^. 



1. Neutralino Decays 



For a neutralino two-body decay to charged slepton and lepton, the flavor-generalized 
matrix element is 



and the decay width is 



(E3) 



(E4) 



22 



For a neutralino two-body decay to sneutrino and neutrino, the matrix element 

* - 



1 , ,9/0 2 ^ ^ 

m,-,. 



and the decay width is 
2. Chargino Decays 

The matrix element for a chargino decay to slepton and neutrino is 
and the decay width is 

T- ( ~- T) ^ Vlrjk^ (2 2 \2 

Xi 

The matrix element for a chargino decay to sneutrino and Icpton is 

= iveM {S^2*PR + S^tPL)v^,{p,), 

and the decay width is 

Xi 

Note that C is the charge conjugation operator. 



^(2) 



3. Slepton Decays 

The slepton decay to neutralino and lepton has a matrix element given by 
M{£i ^ x%) = iueM {PjIIPr + PjfkPr) v^o(pj) , 
and the corresponding decay width is 

+ (-1 - -|o - mi) - 4m,om4 Re {p^P^ 
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For a slepton decay to chargino and neutrino, the matrix element is 
and the decay width is 



2 / \ 2 



TUT ^ ^-Tu,\ = f mj - m"- ) . (E14) 

A slepton decay to hght (neutral) higgs and slepton has a matrix element of 

M{£;^£jh')=ta^f, (E15) 

and the decay width is 

^ ^^^"^ = £!r^''' H' H' -9 ■ (E16) 

The slepton to slepton and Z boson matrix element is 

M{i- ^ i-Z') = iC^f(p,+p.Ye;{pz) , (E17) 

and the decay width is 

Fit ^ tz') = '^^y , A^/^ fm? , m? . m|) . (E18) 

The slepton to sneutrino and W boson matrix element is 

M{i7 ^ i^jW-) = id^ip, + PjYtlipw) , (E19) 

and the decay width is 

|a(1)|2 

Viir ^ = ^^£^^>^"' ■ (E20) 

4. Sneutrino Decays 

The two-body mode for sneutrino decaying to neutralino and neutrino has a matrix 
element of 

Mivi x'j^k) = ioijikUuk(Pk)PRV^o{Pj) , (E21) 

and a decay width 

The sneutrino to chargino and lepton two-body mode has a matrix element of 



Vy+{Pi 



i\ (Pk) {S^IIPr + SI^IPl) (Pj) , (E23) 
24 



and a decay width of 



2 /o-, 2^ 



2 2 

mr, — m 



k 

- <) + 4m-+m,, Re (^j.^^jj))] . (E24) 
The matrix element for sneutrinos decaying to sleptons and W bosons is 

M{v, ^ i-W+) = iC^\, +p.ye;(pw) , (E25) 

and the decay width is 

m - qW^) = TTT^^V^A^/' (m%ml,m\,) . (E26) 

5. Higgs Decays 

Finally, the Higgs decays have trivial matrix elements, so we only present their widths. 
The decay width for heavy Higgs boson to two sleptons is 



while the width to two sneutrinos is 

^ "'""^ = iSk'"' K.">|.0 ■ (E28) 

The decay width for pseudoscalar Higgs to two sleptons is 



and the decay width for charged Higgs to slepton and sneutrino is 

i-^v,) = J^AV^ (ml,,, m|„ mfj . (E30) 



APPENDIX F: THREE BODY DECAYS 

The general three-body decay width is 

body decay TTT^ / ' dE^ dEm\M\ , (Fl) 

where the squared matrix element includes a sum over the outgoing fermion spins. The 
limits are 

= 2(p,-p,)^ [i^i - ^k) (mj - m] + ml + ml,- 2miEk) 
±PkX^{{Pi-Pk)^mlml)] , (F2) 
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and A(x, y, z) = + y'^ + — 2{xy + xz + yz) as in the two-body decays. 

For our purposes, we analytically reduce this to a one- dimensional integral over Ek by 
integrating out Em-, and then we evaluate the integral numerically. We only evaluate the 
three-body modes if there are no kinematically allowed two-body decays^. We could tech- 
nically apply the three-body formulas to the entire slepton sector, but this would in effect 
double count the on-shell decays of heavier sleptons. 

The full form of these slepton three-body modes is too lengthy to include here; complete 
expressions are given in Ref . j22| . At a qualitative level, there are four distinct decay modes: 

r^rf^C r^^^c t^tv^Vm, r^tqmk- (f3) 

The first three decays exist in the flavor-conserving case, occurring through an off-shell 
neutralino or chargino; the addition of flavor violating effects introduces flavor- violating 
couplings and an interference term between the diagrams in the C.~ mode. Flavor 

violating effects also introduce the possibility of radiating an off-shell Z or Higgs, which 
contribute both to the decays to a slepton and a same-generation lepton and anti-lepton 
pair and to a slepton and two quarks. 

APPENDIX G: PROGRAM DETAILS 

Here we outline the structure of the SPICE program. There are three distinct compo- 
nents: the code to impose generalized flavor- violating boundary conditions for SOFTSUSY and 
calculate sleptonic decay widths, an instance of SUSYHIT to calculate non-sleptonic decay 
widths, and a subroutine which merges the information from both files into an output file. 
As previously mentioned, the instance of SUSYHIT is unchanged from its original form. 

1. SOFTSUSY Customization and Sleptonic Decay Widths 

Our customization of the SOFTSUSY code has two components — renormalization bound- 
ary conditions with both gauge and gravity mediated contributions and calculation of decay 
widths for slepton flavor violating decays. 

We have not altered the SOFTSUSY renormalization algorithms; however, since the lepton 
mass matrix is in general non-diagonal, the mixing between lepton gauge and mass eigen- 
states must be rotated out of the lepton sector at the high scale to conform to SOFTSUSY's 
boundary conditions. This is done by first diagonalizing the lepton mass matrix 

^L^rriE^/R, = h^ui'^^m^^U^'/^kR^ = h.m'^jR, , (Gl) 

where mE^^ is the original lepton mass matrix and fn^^, is the diagonalized lepton mass 
matrix. 

The interaction terms in the Lagrangian have the form 

iljL^ = ut^^ul^l^tliL, , (G2) 



We always ignore two-body decays to gravitino. See Sec. |TV]for details. 
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with similar expressions for the right-handed slepton and sneutrinos. Then the lepton mixing 
may be absorbed into the slepton mixing at the high scale, whereby 



Tri^L'f) _^ JT{h)'[TTitL) 

again with similar expressions for right sleptons and sneutrinos. 

Here we take advantage of the negligible left-right mixing terms at the high scale to 
decouple mixing of left and right sleptons, though the procedure does work in general if 
the 3x3 matrix associated with left-right mixing is also rotated in a similar fashion. The 
diagonalized lepton mass matrix and the rotated slepton and sneutrino mass matrices are 
then used as the high-scale boundary conditions for SOFTSUSY's renormalization routine, and 
is found in the gaugegravityBcs method. 



2. Sleptonic Decay Calculations 

The new calculation of flavor violating decays is encompassed in three classes: 
DecayFormat, DecayCalc, and ThreeBodylntegral. Calculations of decay widths is per- 
formed in the DecayCalc class, with the help of the ThreeBodylntegral s class for three- 
body decays. The DecayFormat class handles output to the intermediate file and interfaces 
with the main program. 



a. DecayCalc class 



Method 


Purpose 


Lambda 


Common phase space function 


DecayCalc 


Loads the slepton masses and mixing matrices 




from a passed FlavourMssmSof tsusy object. 


Alpha, Beta[l-2] , 


Methods which return the value of the 


Gamma, Delta [1-2] , 


corresponding coupling coefficient for a given 


Sigma [1-4] , Rho, 


set of mass eigenstates. The methods are in 



Zeta[l-2] direct correspondence to the couplings in the 

Lagrangians presented in App. O 

Calc [ . . . ] Calculates the decay width (in GeV) of the given 

decay for a given set of mass eigenstates. The 
calculation is performed according to the decay 
widths given in App. [E]and App. [Fl 

Output [ . . . ] Decays Methods which loop through all decays of a 

particular particle type and passes decays 
widths to the DecayFormat object for output. 

IntegralCoef f icient Methods to produce the beta coefficients 

CrossCoef f icient which weight the three-body decay integrals. 
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DecayCalc loads the slepton mass and mixing matrices from a passed 
FlavourMssmSof tsusy object, which has aheady undergone the renormahzation rou- 
tine. The slepton and sneutrino masses are contained in six and three component vectors, 
respectively, and the corresponding rotation matrices are 6x6 and 3x3 orthogonal 
matrices. The rest of the class is dedicated to calculations of decay widths. The individual 
decay widths are calculated from the methods named Calc[. . .] , which take in the integer 
index of the mass eigenstates associated with the decaying and daughter particles. Here i is 
the index of the decaying particle, while j, k, and m are the indices of the daughter particles 
in order. Finally the decays are organized in the methods named Output [. . .] Decays, 
which each contain all the decays of the appropriate parent particle with loops over all 
decaying and daughter particle mass indices. The remainder of the class is composed of 
helper functions, including functions corresponding to the functions Lagrangian coefficients 
in App. O Additional helper methods are included in the file DecayMath. cpp, which are 
primarily functions which perform the integration for three-body decays. 



b. ThreeBodylntegrals Class 



Variable 


Description 


start, end 


Upper and lower integration limits calculated 




from the currently saved mass ratios 


r[...] 


Dimensionless ratios of the masses of the 




daughter particles and the mass of the initial 




slepton 


Method 


Purpose 


Lambda 


Common phase space function 


F, Fl, F2 


Common functions used in the integrals 


Simpsons 


Method to perform an integral of a given 




one-dimensional function. The integration 




limits are stored class variables. 


Integral [1-6] 


Six integrals corresponding to the standard 


Crosslntegral [1-8] 


three body decay mode and eight corresponding 




to the intereference term. 


setMasses 


Sets the mass ratios and integration limits. 


evaluatelntegrals 


Methods which group together sets of 


evaluateCrossIntegrals 


integrals and then take a sum using a 


evaluateCharginoIntegral vector of coefficients as weights 



The ThreeBodylntegrals helper class is used to calculate the values of integrals for the 
three body decays. The class takes in particle masses using the function setMasses, which 
sets common integration limits and mass ratios to use for multiple integrals. A summation 
over several integrals (or just one for the chargino) can then be performed using the methods 
evaluatelntegrals, evaluateCrossIntegrals, and evaluateCharginoIntegral. Coeffi- 
cients for the integrals are passed into these methods to properly weight the sum. 
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c. DecayFormat 

Method Purpose 

PDGCode Function which returns the PDG code of a 

particle given its type and mass eigenstate 
OutputParticle Output particle code, mass, and lifetime 

OutputParticleData Output a table of particle masses and lifetimes 

TwoBodyDecay Output the decay width for a two-body decay 

ThreeBodyDecay Output the decay width for a three-body decay 

OutputDecayTable Output a table of decay modes with decay widths 

for decays involving sleptons and sneutrinos 
OutputSusyParams Outputs SUSY parameters to be appended to the 

end of HERWIG input file 
OutputlntermediateFile Outputs the full intermediate file containing 

sleptonic decays 

Finally, the class DecayFormat is devoted to constructing the intermediate file to send 
to the FileCombine part of the code. The DecayFormat class creates an output stream at 
declaration time, and its two methods TwoBodyDecay and ThreeBodyDecay are used by the 
DecayCalc class to output two- and three-body decays, respectively. The remaining methods 
are used to output the particle decay table and HERWIG and SLHA model parameters. 

d. Adding a Decay 

The easiest possible alteration to the program is the addition of another decay mode. To 
add a new decay mode, first add a new function to produce the appropriate decay width in 
corollary to the functions Calc[. . .]; this functions should take in the appropriate indices 
for the mass eigcnstatcs. Then place a set of loops in the appropriate Output [ . . . ] Decays 
function, which calculates the decay width and, if it is non-zero, outputs the decay width 
through either the TwoBodyDecay or ThreeBodyDecay method of the stored DecayFormat 
object in the DecayCalc class. From this, the decay mode will be added to the intermediate 
file and be automatically included in the final HERWIG and SLHA output files. 

3. FileCombine Program 

The FileCombine program merges the intermediate file from the sleptonic decay calcula- 
tion and the SUSYHIT output. The program is composed of three classes: DecayMode stores 
the particle ID's and decay width of one decay mode. Particle stores all the information, 
including a vector of all decay modes, for a single particle. Model contains a map of all 
particles, various model parameters, and both input and output functions. 
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a. DecayMode Class 




Method 


Purpose 


IsSleptonic 


Returns true if the parent or any daughter 
particle is a slepton or sneutrino 


addDaughter 


Add another daughter particle 


OutputHerwigDecayMode 

OutputSLHADecayMode 

0utputSLHA2DecayMode 


Output the decay mode in HERWIG format 

Output the decay mode in SUSY Les Houches format 

Output the decay mode in SUSY Les Houches 2 format 



The DecayMode class stores all information for one decay. This includes the parent particle 
HERWIG ID, the number of daughter particles, the PDG and HERWIG ID's of the daughter 
particles, and the decay width of the decay mode (in GeV). The class also contains methods 
which output the mode in HERWIG and SLHA(2) format. 



b. Particle Class 



Method 


Purpose 


AddDecayMode 


Add given decay mode to the decays for the particle 


FixWidth 


Determine the total decay width based on all currently 
stored decay modes and produce branching fractions 



OutputHerwigDecays Outputs the decay table in HERWIG format 

Output SLHADe cays Outputs the decay table in both SLHA and SLHA2 format 

The Particle class contains the mass, total decay width, and decays for a given particle. 
The decays are stored as a vector of DecayMode objects. 



c. Model Class 



Variable 


Description 


particleMap 


Map object containing stored particle objects. The map 




is indexed by HERWIG Id's 


Method 


Description 


PDGToHerwig 


Converts from PDG to HERWIG particle code 


OnShell 


Checks whether a passed decay mode is on shell 


Readlntermediate 


Read in the model and decay information from the 




intermediate file produced by SPICE 


ReadLesHouchesFile 


Read in decays from the Les Houches file produced 




by SUSYHIT 



OutputHerwigFile Produces a HERWIG input file from the 

stored model parameters and decays 
OutputLesHouchesFile Produces a SLHA and SLHA2 output file. 



This method copies the previous decay file up 
the decays, then produces a new decay table. 
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The Model class contains a map of particle objects, stores the generic model parameters, 
and handles the input/output functions of FileCombine. The Readlntermediate function 
reads in particle mass information, decays, and model parameters from the intermediate 
file; ReadLesHouchesFile fills the rest of the decay table from the SUSYHIT output file, 
discarding off-shell decays. The Combine function is called from outside the class to force 
a consistent order (since the particle mass information is read from the intermediate file). 
The two output functions produce HERWIG input and SLHA and SLHA2 format files. 
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